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Abstract 

This paper shows the application of Genetic 
Algorithm to Reentry trajectory optimization 
problem in the presence of constraints such as 
dynamic pressure and heat flux. The coupled non 
linear equations of motion of reentry body are 
difficult to solve using conventional methods like 
finding costate equations and then applying 
standard optimization methods. As an alternate, 
stochastic way of optimization through genetic 
algorithm is proposed. The results for a typical 
Reusable Launch Vehicle - Technology 
Demonstrator are in good match with that of Space 
Shuttle. 


Introduction 

The concept of Reusable Launch Vehicle 
(RLV) is evolved in order to minimize cost per unit 
payload. Nearly 70% of cost of a launch vehicle 
lies in structure and avionics system, recovering the 
same and using it for further flights (after 
refurbishing) will be less expensive. The 
technologies like TPS material, structure and 
avionics system for a full fledged RLV are highly 
challenging. As a cost effective measure, a smaller 
version of RLV is proposed called as Technology 
Demonstrator (RLV-TD) as given in Fig.1 which 
will validate the above mentioned technologies in a 
step by step manner. 


Nmenclature 
Cp 7 drag coefficient 
Cpo = constant 
CL = lift Coefficient 
Cho = C,_ as a function of a, /rad 
D = drag force, N 
Fr — thrust force, N 
Zo 7 gravity at earth surface, m/s” 
g = gravity at some altitude, m/s” 
h = altitude, m 
Ty = specific impulse, N s/kg Fig.1 RLV — Technology Demonstrator 
k = constant 
L = lift force, N In the reentry trajectory design one has to 
LH = local horizontal plane ensure that with minimum fuel consumption the 
M = total mass, kg desired landing site is reached satisfying 
Q = dynamic pressure, Pa constraints such as dynamic pressure and heat flux. 
q = heat flux, W/cm? Traditionally reentry trajectory optimization is 
To = radius of earth, m carried out by non dimensionalising the equations 
r = radial vector, m of motion! and then solving it through gradient 
Ry = nose radius, m method using algorithms such as SGRA” 
S = area, nm Brysons*. As the equations of motion of a reentry 
U = control vector body are coupled non linear equation, it gives 
V = velocity, m/s complicated co state equations which are difficult 
a angle of attack, deg to solve and simulate. As an alternate, evolutionary 
B ea eoactint approach is proposed for this type OF problems. 
= : Genetic Algorithms(GA)” are — search 

: Ee ce ook deg algorithms based on the mechanism of natural 

, ; 3 selection. They lay on one of the most important 
Po density at earth surface, kg/m par oe \ ‘ 

: : 3 principles of Darwin : survival of the fittest. 
, 7 density at some altitude,kg/m Globally we make use of the population, submitted 
8 = longitude, deg to many transformations. After some generations, 
S = bank angle, deg the population endues no more; the best individual 
Vv = heading angle, deg represents the optimal solution. GAs do not require 
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derivative or continuity of a function to be 
evaluated. They are simple, fast converging and 
easy to implement and simulate. 


Equations of motion of a reentry body 


The equations of motion for a lifting trajectory are 
given’” by (Fig.2) 
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Fig.2 Definition of various parameters 


The control vector is U=[a,o]. The drag, the lift, 
the density and the gravity acceleration are 
assumed to have standard form : 
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The density is assumed to be zero above 120km 
altitude. The heat flux is ae by 


Genetic Algorithm 


They are different from more normal optimization 
and search procedures in following ways”” 

- GAs work with coding of the parameter set, not 
the parameter themselves 

- GAs search from a population of points, not a 
single point. They climb many peak in parallel and 
not go point by point. 

- GAs use objective function information, not 
derivatives 

- GAs use probabilistic transition rules, not 
deterministic rules. 

Traditionally in GA, the number of 
variables in the objective function are coded into 
binary strings (0 and 1). Each bit is a gene and the 
complete string is called a chromosome. An initial 
population (of chromosome) is chosen randomly. 
Each chromosome represents a potential solution to 
the problem. At each generation, the chromosome 
of the population are evaluated from their fitness 
value and reproduced further using mutation and 
crossover. 

In the present study, each chromosome 
vector is represented with floating point numbers. 
It has been shown _ that® floating point 
representation is faster, more consistent from run to 
run and provided a_ higher precision. Its 
performance can be enhanced by special operators 
(like boundary mutation) to achieve high accuracy. 


Problem Statement 


Minimize, the control force 
J=J]U*(t) dt 

and meet constraints on Q and q. 
The integration of equations of motion is carried 
out at 1 s interval upto 300 s. The problem is 
discretized into 300 steps [t),ts,..ts3o00] where we 
have to determine control vector U[a,o] so as to 
satisfy constraints. The initial values of state 
vectors are specified. In the present study, cross 
range optimization is not carried out and as such o 
= 0 is assumed. Also all the propellant is utilized 
from t=0 to t=200s. 


Input Data 


Mass at orbit = 1750 kg 

Mass of propellant = 300 kg 

Isp = 2696.8 N s/kg, S = 0.785 m’, Ry = 0.25m 
Po = 1.23 kg/m’, B = 1/7320, r, = 6378000m 
g, = 9.80665 m/s” 

Initial [h, V,y,0,o,w] is [150000, 1500,0,0,0,0] 
The aerodynamic model is assumed to be 
Cp = Cp, + kc, 

Cp, = 0.4, k = 0.18 

Cc. = Cia -Q 

Cia = 7.0 /rad 
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Results and Discussion 


The constraint boundaries for RLV-TD is 
given in Fig.3. The boundary where lift force is 
balancing gravity minus centrifugal force is called 
Equilibrium Glide Boundary. The vehicle will be 
following heat flux boundary (q) till about 25 km 
thereafter it will go along dynamic pressure (Q) 
boundary . 


Heat flux constraint (< 18.5 W/em’) 


Q constraint (< 45 kPa) 
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Fig.3 Boundary constraints 


The unconstrained Q is > 300 kPa for 
RLV-TD. Using o as controlling parameter, Q can 
be controlled to Q < 50 kPa and Q < 45 kPa as 
given in Fig.4. The control angle history is given in 
Fig.5 for the two cases. 
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Fig.4 Q constraint only 


In the next set of optimization two set of 
constraints were specified 

(1) Q<45 kPa 

(2) q < 18.5 W/cm? 

The altitude, velocity, flight path angle 
and heat flux plot are given in Fig.6, Fig.7, Fig.8 
and Fig.9 respectively. The control angle history is 
given in Fig.10. It is touching maximum 
permissible value of 30deg at 170 s. 
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Fig.6 Altitude history for both Q & q constraint 
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Fig.7 Velocity history for both Q & q constraint 
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Fig.9 Heat flux history history 


30 


20 


Control Angle a (deg) 


10 


100 200 300 


Time s 


Fig.10 Control history for both Q & q constraint 


For the reference trajectory, control U can 
be given as a function of time (Similar to Space 
Shuttle*)as shown in Fig.11. The trend matches 
well with the reference trajectory of Space Shuttle’®. 
The drag acceleration profile as given in Fig.12 
also matches well (in trend) with that of space 
shuttle. 


Summary 


Reentry trajectory optimization problem is solved 
using genetic algorithm for a typical RLV- 
Technology Demonstrator. Dynamic pressure and 
heat flux are used as constraints. The results are 
comparable to that of Space Shuttle. 
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Fig.11 Reference Control 
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Fig.12 Drag acceleration 
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